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A large  data  base  is  collected  from  the  open  literature  and  private 
sources  concerning  key  side  force  characteristics.  Linear  regres- 
sion techniques  are  used  to  establish  predictive  equations  which  can 
be  used  by  the  designer. 
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INTRODUCTION 


Bodies  with  pitch  plane  symmetry  experience  steady  or  quasi-steady 
side  loads  at  hiqh  angles  of  attack  in  the  subsonic- transonic  range. 
These  forces  generally  occur  between  20  to  60  degrees  incidence  and 
can  be  traced  to  the  formation  of  an  asymmetric  vortex  pattern  in  the 
leeward  flow  tield.  The  magnitude  of  the  side  force  has  been  observed 
to  exceed  that  of  the  norma]  force  on  short  bodies  (1^  < 4).^ 


To  date,  side  forces  have  been  the  subject  of  several  studies, 
which  have  determined  that  the  magnitude  of  these  loads  decreases  with 
increasing  Mach  number,  decreasing  noc>e  fineness,  and  the  introduction 
of  bluntness.  No  conclusions  have  been  offered  concerning  the  influ- 
ence of  Reynolds  number,  and  the  addition  of  grit  has  been  observed  to 
increase  these  forces  in  some  cases,  and  decrease  it  in  others. 


A quantitative  description  of  side  forces  has  not  emerged,  these 
loads  being  extremely  sensitive  to  factors  that  are  difficult  to  con- 
trol, such  as  extremely  small  irregularities  in  the  model  geometry 
and  seemingly  insignificant  changes  in  test  conditions  and/or  pro- 
cedures. Rolling  an  axis-symmetry  model  drastically  alters  the 
measured  side  force  and  the  accompanying  vortex  pattern  as  shown  in 
Figures  1 and  2.  The  pictures  of  Figure  1 are  unpublished  results 
obtained  at  NAVSURFKPNCEN  while  the  data  of  Figure  2 is  from 
reference  5.  The  cause  of  this  change  has  been  traced  by  some 
observers  to  irregularities  in  model  construction  on  the  order  of 


1.  Keener,  Earl  R.  and  Chapman,  G.  T. , "Onset  of  Aerodynamic  Side 
Forces  at  Zero  Siceslip  on  Symmetric  Forebodies  at  High  Angles 
of  Attack,"  AIAA  taper  74-770,  1974. 

2.  Pick,  G.  , "Investigation  of  Side  Forces  on  Ogive-Cylinder  Bodies 
at  High  Angles  of  Attack  in  the  M = 0.5  to  1.1  Range,"  AIAA  Paper 
No.  71-570,  Jun  1971. 

3.  Fleeman,  E.  L.  and  Nelson,  R.  C.,  "Aerodynamic  Forces  and  Moments 
on  a Slender  Body  v;ith  a Jet  Plume  for  Angles  of  Attack  up  to  180 
Degrees,"  AIAA  Paper  No.  74-110,  1974. 

4.  Clark,  W.  H. , Peoples,  J.  R.  and  Briggs,  M.  M.,  "Occurrence  and 
Inhibition  of  Large  Yawing  Moments  During  High  Incidence  Flight 
of  Slender  Missile  Configurations,"  J.  Spacecraft,  10,  8,  1973, 
pp.  510-519. 

5.  VJardlaw,  A.  B.,  Jr.,  ''The  Prediction  of  Normal  Force,  Pitching 
Moment  a.nd  Yaw  Force  on  Bodies  of  Revolution  at  Angles  of  Attack 
Up  to  50  Degrees  Using  a Concentrated  Vortex  Flaw-Field  Model", 
NOLTR  73-209,  Oct  1973. 
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FIG  1 VARIATION  OF  THE  LEEWARD  FLOWFIELD  WITH  CHANGE  IN  THE  MODEL  ROLL  ORIENTATION. 
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DATA  TAKEN  FROM  REFERENCE  5.  (M  = 0.9) 
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FIG.  2 Cy  AS  A FUNCTION  OF  a AT  SEVERAL  DIFFERENT  ROLL  ORIENTATIONS. 


NS'WC/WOL/TR  75-176 


m?3hining  tolerance; ® however,  the  fact  that  side  force  measurements 
taken  at  roll  angles  differing  by  180  degrees  are  not  necessarily 
mirror  images  of  one  another  (see  Figure  2)  indicates  that  other  fac- 
tors are  important.  In  particular  it  raises  the  questions  as  to 
whether  side  force  magnitude  is  facility  dependent. 


The  definition  of  the  high  angle  of  attack  flowfield  is  further 
complicated  by  wake  instabilities  which  have  been  cbserved  to  be  inter- 
spersed with  the  steady  asymmetric  mode.^  Although  the  steady  mode 
appears  to  be  the  dominant  one,  a completely  steady  asymmetric  flow 
field  and  aerodynamic  load  model  is  an  idealization  of  the  situation. 
The  results  of  reference  6,  which  at  present  stand  alone,  indicate  a 
predominant  unsteady  mode  at  angles  of  attack  as  low  as  10  degrees. 


In  the  current  report  a quantitative  description  of  key  side  force 
characteristics  is  developed  by  applying  linear  regression  techniques 
to  a collected  data  base.  The  data  base  is  constructed  from  informa- 
tion available  in  the  open  literature  and  from  private  sources  con- 
cerning the  key  features  of  maximum  observable  side  force,  angle  of 
attack  at  which  it  occurs,  and  the  lowest  angle  of  attack  at  which 
a side  force  is  observed.  As  a prelude  to  applying  linear  regression 
methods  it  is  necessary  to  postulate  a set  of  independent  variables 
that  is  thought  to  influence  each  of  these  quantities.  This  could  be 
approached  from  the  dimension  analysis  point  of  view;  however,  t-his 
produces  too  many  possible  independent  variables.  Instead,  an  approx- 
imate model  based  on  the  crossflow  analogy  will  be  constructed  to 
predict  side  force  and  its  onset  angle.  The  results  of  this  analysis 
in  conjunction  with  a detailed  study  of  the  data  base  will  be  used  to 
postulate  a set  of  independent  variables  related  to  each  of  the  above 
quantities.  The  linear  regression  procedure  will  determine  a least 
square  fit  for  each  of  these  quantities  using  only  those  independent 
variables  which  are  statistically  significant. 


PREDICTIVE  MODEL  BASED  ON  THE  CROSSFLOW  ANALOGY 


The  crossflow  analogy  assumes  that  forces  on  a body  at  angle  of 
attack  are  related  to  those  acting  on  a cylinder  in  crossflow  and  has 
been  extensively  applied  to  the  prediction  of  normal  force.  In  its 
simplest  form  the  steady  state  cylinder  drag  corresponding  to  and 
Re^  is  used  in  computing  this  quantity.®'^  A more  accurate  approach 


6.  Gowen,  F.  E.  and  Perkins,  E.  W. , "Study  of  the  Effects  of  Body 
Shape  on  the  Vortex  Wakes  of  Inclined  Bodies  at  M=2,"  RM  A53117, 
1953,  NACA. 

7.  Thomson,  K.  D.  and  Morrison,  D.  F. , "The  Spacing,  Position  and 
Strength  of  Vortices  in  the  Wake  of  Slender  Cylindrical  Bodies 
at  Large  Incidence,"  J.  Fluid  Mechanics,  50  4,  1971,  pp.  751-783. 

8.  Allen,  H.  J. , "Estimation  of  Forces  and  Moments  on  Inclined 
Bodies  of  Revolution  of  High  Fineness  Ratio,"  NACA  RM  A9I26, 

1949. 

9.  Jorgensen,  L.  H.,  "Prediction  of  Static  Aerodynamic  Characteris- 
tics for  Space- Shuttle-Like  and  Other  Bodies  at  Angles  of  Attack 
from  0°  to  180°,"  NASA  TN  D-6996,  Jan  1973. 
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assumes  that  the  crossflow  plane  is  swept  down  the  length  of  the  inci- 
dent body  at  the  uniform  rate  U cos  a.  Distance  along  the  body  is 
related  to  time  in  the  case  of  a cylinder  by: 


t = x/(Ucosa) 


(1) 


At  any  point  along  the  model  length  the  steady  state  drag  value  is 
replaced  by  that  of  an  impulsively  started  cylinder  at  its  equivalent 
time  t.  This  approach  was  originally  suggested  by  Kelly^O  and  has 
since  been  used  in  various  forms  by  others. 

In  considering  a predictive  model  for  side  force,  it  is  necessary 
to  use  the  second  approach  since  the  lift  on  a cylinder  is  unsteady 
while  Cy  on  a body  at  incidence  is  steady.  Hence: 

L 

Fy  = J'  (t=x/Ucosa,  dx  (2) 

0 


Ideally  Fl  would  be  based  on  data  for  an  impulsively  started, 
expanding  cylinder.  Unfortunately,  this  type  of  information  is  not 
available  so  time  independent  root  mean  square  values  are  used. 

Substitution  of  the  unsteady  lift  of  a cylinder  in  crossflow 
into  the  right  hand  side  of  Equation  (2)  results  in  a force  distri- 
bution which  varies  periodically  from  positive  to  negative  along  the 
length  of  the  body.  Since  the  frequency  for  the  shedding  of  consec- 
utive vortices  on  a cylinder  is: 


f = 2US/D 


(3) 


the  distance  along  the  body  sustaining  a force  of  constant  sign  is: 

X = D/(2  S tana)  or  1 = 1/(2  S tana)  (4) 

s s 


Hence,  the  maximum  unbalanced  side  force  is  obtained  by  integrating 
Equation  (2)  between  the  limits  C and  C + Xg  where  C mar)cs  the  dis- 
tance in  X at  which  the  side  force  changes  sign.  Such  an  interval 
containing  a force  of  constant  sign  will  be  referred  to  as  a cell. 
In  terms  of  yaw  and  lift  force  coefficients  Equation  (2)  becomes: 


10.  Kelly,  H.  R. , "The  Estimation  of  Normal  Force,  Drag,  and  Pitching- 
Moment  Coefficients  for  Blunt  Based  Bodies  of  Revolution  at  I-arge 
Angle  of  Attac)c,"  Journal  of  Aeronautical  Sciences,  Vol.  21, 

No.  8,  pp.  549-555. 
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4D  sin  a 


/ 


4- 


D/(2  S tana) 


(5) 


Here  D represents  an  average  value  for  D over  the  length  of  the  shed- 
ding interval.  is  generally  expressed  by  its  root  mean  square 
value  and  assuming  that  this  measurement  is  generated  either  by  sine 
or  square  wave  results  in: 


.3  sin  (2  a)  C^. 

rms 

S 


(6) 


The  extent  of  the  available  experimental  values  for  and  S are 
given  in  Figures  3 and  4 respectively  and  have  been  taken  from 
references  11  through  14.  The  shaded  area  represents  the  range 
in  measured  values  for  given  test  conditions.  Although  the  incompres 
ible  case  is  well  documented,  very  little  information  is  ai'ailable 
in  the  transonic  or  supersonic  regime.  Transonic  measurements,  taken 
at  supercritical  Reynolds  numbers,  indicate  that  as  the  Mach  number 
is  increased  past  the  critical  point,  Cj^  changes  noticeably  assum- 
ing a new  value  that  is  similar  to  that  for  the  laminar  case.^^  This 
same  type  of  behavior  is  characteristic  of  the  cylinder  drag.  For 
supersonic  conditions  no  direct  measurements  are  available;  however, 
Schlieren  photographs  show  a symmetric  wake  structure. Based  on 
this  evidence  it  is  likely  that  Cj,  = 0 under  these  conditions.  In 
Figure  4,  the  measurements  of  reference  7 for  S on  a body  at  incidence 
with  laminar  boundary  layer  are  shown.  These  agree  well  with  those 
taken  on  a cylinder  in  crossflow. 

The  preceding  model  does  not  provide  a prediction  of  onset  angle. 
To  develop  such  an  expression  the  sequence  of  events  accompanying  the 
formation  of  the  flow  field  about  an  impulsively  started  body  is 


Hi" Jones,  G.  W. , Cincotta,  J.  J. , and  Walker,  R.  W. , "Aerodynamic 
Forces  on  a Stationary  and  Oscillating  Circular  Cylinder  at  High 
Reynolds  Number,"  Oct  1968,  NASA  TR  R-300. 

12.  Bearman,  P.  W. , "On  Vortex  Shedding  From  a Circular  Cylinder  in 
the  Critical  Reynolds  Number  Regime,"  J.  Fluid  Mech.,  Vol.  37, 

Part  3,  pp.  577-585,  1969. 

13.  Richter,  Andreas,  Flow  Force  on  Rigid  Circular  Cylinder  Between 
Parallel  Walls,  (in  German),  PH.D.  Dissertation,  Karlsruhe,  1973. 

14.  Gowen,  F.  E.  and  Perkins,  E.  W.,  "Drag  of  Circular  Cylinders  for  a 
Wide  Range  of  Reynolds  Numbers  and  Mach  Num±)ers,"  NACA  RM  A52C20, 
Jun  1952. 
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considered.  Sarpkaya's  data^^  taken  in  incompressible  flow  over 
the  Reynolds  number  range  1.5  x 10^  to  1.2  x 10^  indicates  that  the 
flowfield  remains  symmetric  for  an  initial  period  t*g  = Ut/D  = 4 , 
Generalizing  to  a body  at  incidence  and  using  D in  place  of  D to 
account  for  the  changing  body  diameter  near  the  nose; 


t*  = X (tan  a )/D  = 4 
s o 


and  (7) 

= tan“^  (4/1^) 


Equation  (7)  does  not  indicate  the  effect  of  Mach  number  on  uq. 

Based  on  slender  body  crossflow  theory  it  might  be  expected  that  this 
effect  would  not  be  of  great  importance  since  Oq  is  often  relatively 
small. 

To  produce  a side  force  estimate  continuous  in  a,  it  is  necessary 
to  define  ap  as  the  lowest  angle  of  attack  at  which  a complete  cell 
exists  on  the  body.  To  determine  ap,  It  “ Is  is  used  in  place 
of  X in  Equation  (7); 


If  a is  between  Uq  and  ap,  only  the  following  fraction  of  the 
cell  is  present: 

-—  = 2S  [^1^  tana  -4j 

If  a sine  force  distribution  is  assumed  along  the  length  of  the  body 
at  incidence,  the  value  of  Cym  calculated  in  Equation  (6)  must  be 
multiplied  by  the  correction  factor  F = (1/2  - 1/2  cos  [Trl/lg])  in 
the  angle  of  attack  range  (Xq  to  ap. 

The  predictions  of  the  crossflow  model  are  summarized  in  Figures  5 
and  6.  In  these  graphs  it  has  been  assumed  that  the  maximum  incidence 
at  which  a steady  side  force  occurs  is  60  degrees.  The  value  of  Cl 
used  is  indicated  by  the  dotted  line  in  Figure  4 and  3 is  taken  as 
.17  .2.  The  calculated  curves  in  Figures  5 and  6 reflect  the  fact 


Xm  Sarpkaya , T.,  "Separated  Flow  About  Lifting  Bodies  and  Impulsive 
Flow  About  Cylinders,"  i lAA  Journal,  Vol.  4,  No.  3,  1966, 
pp.  414-420. 
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that  the  maximum  side  loads  occur  at  a = 45**  for  M < .46  as  su«;;- 
gested  by  Equation  (6) . At  higher  Mach  numbers,  variations  in  Cj, 
dominate  over  the  angular  dependence  of  this  equation  and  Cym  is 
located  at  the  angle  of  attack  at  which  Mg  - .55  and  Cj,  has'^its  peak 
value.  In  the  incidence  range  olq  to  ap,  F decreases  rap_idly  with 
decreasing  a and  the  peak  side  force  occurs  at  ap.  For  1^  < 8, 

Up  is  greater  than  45**  which  explains  the  unique  characteristics  of 
the  = 4 curve. 

The  crossflow  model  indicates  that  Cvm  is  a function  of  M and  Re 
for  < Merit-  The  model  length,  1^/  is  only  important  at  high 
Mach  numbers  and  on  very  short  models,  while  nose  geometry  is 
significant  only  for  models  without_afterbodies.  In  this  case  Cy^ 
decreases  with  decreasing  Ij^  since  D becomes  smaller  under  these"^ 
conditions.  The  variable  am  is  dep  .ndent  on  M and  It  while  Uq 
is  a function  of  It  only. 


DATA  BASE 

The  collected  data  base  is  given  in  Table  1 and  is  directly  read 
from  references  1 through  3,  5,  and  16  through  24.  The  tabulated 
quantities  Cyjp,  am/  and  a^  represent  the  absolute  value  of 
the  maximum  side  force,  the  angle  of  attack  at  which  it  occurs,  and 
the  lowest  angle  at  which  a side  force  is  measured,  respectively. 

IS'."' Smi'th,  L.  H.  and  Nunn  R.  H., , "Aerodynamic  Characteristics  of  an 
Axisymetric  Body  Undergoing  a Uniform  Pitching  Motion,"  Naval 
Postgraduate  School,  NPS-59Nn75021,  Feb  1975. 

17.  Coe,  P.  L.,  Jr.  and  Chambers,  J.  R. , "Asymmetric  Lateral-Direc- 
tional Characteristics  of  Pointed  Bodies  of  Revolution  at  High 
Angles  of  Attack,"  NASA  TN  D-7095,  Nov  1972. 

18.  Jorgensen,  L.  H.  and  Nelson,  E.  R. , "Experimental  Aerodynamics 
Characteristics  for  a Cylindrical  Body  of  Revolution  with  Various 
Noses  at  Angles  of  Attack  from  0 to  58“  and  Mach  Numbers  from  .6 
to  2.0,"  NASA  TMX-3128,  Dec  1974. 

19.  Jorgensen,  L.  H.  and  Nelson,  E.  R. , "Experimental  Aerody ramie 
Characteristics  for  Bodies  of  Elliptic  Cross  Section  at  Angles 
of  Attack  from  0 to  58“  and  Mach  Numbers  from  .6  to  2.0," 

NASA  TMX-3129,  Feb  1975. 

20.  Daniels,  P.,  "Minimization  of  Lock-In  Roll  Moment  on  Missiles 
via  Slots,"  NSWe/DL  Technical  Report  TR-3250,  Jan  1975. 

21.  Fellows,  K.  A.,  "Tests  on  an  Ogive  Cylinder  up  to  Very  High 
Incidences  at  Transonic  Speeds,"  Aircraft  Research  Association 
Limited,  Bedford,  Great  Britain,  Model  Test  Note  M.61/1. 

22.  Brown,  R.  C.,  "On  the  Asymmetrical  Aerodynamic  Forces 
of  Slender  Bodies  of  Revolution,"  Proceedings  of  BOWACA 
Meeting,  1965. 

23.  "Unpublished  Tests  on  Cone-Cylinder  at  Incidences  up  to  70," 
provided  by  K.  D.  Thomson,  Australian  Research  Laboratory, 

Nov  1971. 

24.  Fidler,  J. , Unpublished  Data  Taken  in  NSRDC  Transonic 
Tunnel,  Private  Communication,  Mar  1973. 
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In  cases  where  two  side  force  peaks  of  nearly  equal  magnitude  are 
present r the  second  peak  and  the  angle  of  attack  at  which  it  occurs 
are  noted  and  included  under  the  headings  of  C'yjj,  and  The 

onset  angle  is  taken  to  be  the  largest  angle  or  attack  at  which  zero 
side  force  is  measured.  Not  all  the  entries  in  Table  1 have  values 
for  Co  and  csjji-  In  some  cases  these  values  could  not  be  read 
clearly,  and  in  others  that  data  was  taken  from  charts  showing  Cyj^j 
only. 

Figures  7 and  8 indicate  that  M or  M^,  is  an  important  determinant 
of  Cym  particularly  in  the  high  subsonic  regime  where  the  crossflow 
Mach  number  is  supercritical.  The  general  trend  is  for  a decreasing 
value  of  Cym  with  increasing  Mach  number.  Smaller  Ij^  values 
result  in  decreasing  values  of  Cyu,  as  is  illustrated  in  Figure  9. 

A detailed  examination  of  the  data  base  indicates  that  this  variation 
is  most  pronounced  at  high  subsonic  Mach  numbers.  The  relation 
between  Re  and  Cy^i  is  explored  in  Figure  10  which  suggests  a sma_l 
decrease  at  higher  Reynolds  numbers. 

The  differences  in  the  relations  governing  blunt  and  sharp  tan- 
gent ogive  of  bodies  are  seen  by  examining  Figures  7 through  9. 
Dependencies  are  much  more  clearly  defined  for  sharp  bodies,  particu- 
larly for  1^.  The  relation  between  Cyj^  and  rj^  on  blunt  bodies  is 
shown  in  Figure  11  and  does  not  indicate  a clear  trend  with  varia- 
tions in  nose  bluntness. 

An  examination  of  Figure  7 shows  that  significant  side  forces 
disappear  for  M > 1.4,  the  only  exception  being  a measurement  taken 
on  a cone.  Another  example  is  available  in  the  literature  showing  a 
small  side  force  on  a cone  at  M = 2.^^  This  suggests  that  cones 
experience  side  force  to  a Mach  number  higher  than  do  tangent  ogives. 
In  Figure  8 it  is  shown  that  side  force  disappears  on  all  classes  of 
bodies  for  M^  > .8.  Thomson  found  that  at  ' .7  the  structure 

of  the  leeward  flow  field  changes  noticeably  with  vortices  forming 
only  near  the  body  nose  and  at  higher  crossflo.?  numbers  disappearing 
altogether. ^ 

Examination  of  the  data  base  ind’ cates  that  am  is  related  to  M 
and  lip  while  is  a function  of  1^.  alone.  Data  for  these  three 

relations  are  presented  in  Figures  12  through  14,  respectively.  In 
all  three  figures  only  values  of  and  Oq  corresponding  to  cases 
where  Cym  i .75  are  included.  An  excessive  amount  of  scatter  was 
found  in  the  remaining  points  due  to  a difficulty  in  reading  them. 


25.  Orlik-Riickerman,  K.  J.,  La  Serge,  J.  C.,  Iyengar,  S.,  "Half 

and  Full  Model  Experiments  on  Slender  Cones  at  Angle  of  Attack," 
Journal  of  Spacecraft  and  Rockets,  Vol.  11,  No.  9,  pp.  575-580, 
Sep  1973. 
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FIG.  8 C.,„  AS  A FUNCTION  OF  = M sin  a „ FOR  ALL  POINTS  IN  THE  DATA  BASE  FOR  WHICH  a„  IS  KNOWN. 
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FIG.  9 C „ AS  A FUNCTION  OF  £„  FOR  ALL  POINTS  IN  THE  DATA  BASE. 
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FIG.  1 1 Cy^  AS  A FUNCTION  OF  rj,  FOFt  ALL  BLUNT  MODELS  (fj,  > 0.005). 
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FIG  13  a„,  AS,  FUNCTION  M FOR  ALL  POINTS  IN  THE  DATA  BASE  WHERE  IS  KNOWN  AND  C„„  >.75. 
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CREATION  OF  A REGRESSION  MODEL 

The  predictive  model  is  to  be  developed  using  Efroymson's  multiple 
regression  analysis  technique^^.  Once  a model  has  been  postulated 
this  method  adds  terms  one  by  one  to  the  regression  equation  taking 
the  term  with  the  highest  correlation  coefficient  first.  Before  this 
variable  is  actually  included  in  the  model  the  F test^^  is  applied  to 
its  variance  contribution  to  determine  whether  it  is  statistically 
significant  at  a specified  F level.  At  the  end  of  each  step,  prior 
to  the  addition  of  a new  variable,  all  variables  in  the  model  are 
checked  to  determine  if  they  are  statistically  significant  and  those 
which  are  not  are  dropped. 

In  order  to  successfully  determine  a predictive  model  using  the 
above  algorithm,  the  correct  variables  must  be  included  in  the  pro- 
posed set.  Although  statistically  insignificant  variables  are  not 
incorporated,  it  is  possible  for  incorrect  terms  to  be  included  as  a 
result  of  peculiarities  in  the  data  base.  Accordingly,  in  order  to 
develop  a predictive  equation  which  is  physically  reasonable,  it  is 
necessary  to  postulate  an  initial  model  which  is  as  accurate  as 
possible. 

A model  of  Gym  is  considered  first.  Consistent  with  predictions 
of  the  crossflow  model  and  the  general  trends  evident  in  the  data  base, 
the  following  model  is  assumed; 

n 

^7  ""b  ^8  (9) 

■^t 

Here  f(Re)  and  g(rb)  are  defined  as  follows; 


^6.  Efroymson,  M.  A.,  "Multiple  Regression  Analysis,  Mathematical 
Methods  for  Digital  Computers,"  Ralston  and  Wilt,  John  Wiley  & 
Sons,  New  York,  1962,  pp.  1^1-203. 

27.  Draper,  N.  R.  and  Smith,  H.,  "Applied  Regression  Analysis," 
John  Wiley  Sons,  1966,  pp.  24-26. 
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f (Re) 


q(^) 


0 Re  < 6(10^) 

1 Re  > 6(10^) 


M-.57 


M > .57 


0 M < .57 


or  M > .57 
and  M < .57 


(10) 


The  form  for  f(Re)  and  q(— ) assumed  is  consistent  with  the  type  of 

relations  suggested  by  the  crossflow  model  where  Re  and  play  a very 
significant  role.  The  variable  1^.  has  been  included  in  inverted  form 
to  insure  that  in  the  limit  In  -»■  “ a finite  force  is  obtained.  The 
parameter  rj^  is  included  since  blunting  is  observed  to  have  a large 
effect  on  Cyjp. 

The  data  base  examination  carried  out  in  the  preceding  section 
establishes  that  the  type  of  body  has  a strong  influence  on  side 
force.  This  suggests  that  a separate  model  be  developed  for  each 
class  of  body,  a procedure  which  is  practical  in  terms  of  the  avail- 
able data  base  only  for  sharp  and  blunt  tangent  ogives.  The  proposed 
model  of  Equation  (7)  is  applied  to  both  types  of  bodies  with  b7  = 0 
for  sharp  tangent  ogives.  Requiring  a 95%  confidence  level^^  for  the 
entering  or  leaving  of  a particular  variable  results  in  the  predictive 
equation: 

(a)  Sharp  tangent  ogives  (rj^  < .005) 


C = 4.20  - 12.5  h(^)  - 1.06  f(Re) 


Standard  Error  of  Cy^  = .65  where: 


h'(^)  = 
n 


. 3 for 


M 


1 - 
n 


> .3 


— T~ 

n n 


< .3 


(11) 
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(b)  Blunt  tangent  ogives  (r^  > .005) 

C = 1.17  - .543  f(Re)  - 13.91  h”  (--) 
ym  -n 

Standard  Error  of  C„_  = .380 

ym 


0 M < 

.57 

M-.57 

M > 

^n 

(12) 


Defining  special  limits  on  M/ln  as  described  by  the  functions  h' 
and  h"  were  ‘juggested  by  an  examination  of  the  residuals. 

As  an  alternative  to  Equation  (11)  a model  is  proposed  using  Me 
and  Rec  in  place  of  M and  Re  respectively  for  tangent  ogives.  Here 
again  a solution  in  the  form  of  Equation  (7)  is  assumed  with  b7  = 0 
and  f(Re)  defined  as: 


f (Re) 


0 Re  < 4,25  (10^)  or  > .5 

c c 

1 Re  >4.25  (10^)  and  < .5 

c c 


This  yields: 


C = 3.65  - 1.019  - 8.39  h (t^) 

ym  cl 


.889  f(Re) 


Standard  Error  of  C, ^ 

ym 

”c 

M 1 C - 

h(-,r2)  " 


= ,70 


n 


.3 

M 


n 


M 

n 


(13) 


Although  this  equation  does  not  provide  as  good  a fit  as  Equation  (9) , 
it  can  be  used  to  generate  Cyjj,  as  a function  of  a. 
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The  variable  am  is  related  to  M and  It  by  both  the  crossflow 

model  and  the  data  base.  Accordingly  the  following  model  is  postulated 

for  a^s 
m 


a_ 


m 


+ bj^M  + b2M 


2 


(14) 


Here  again  It  is  included  in  inverted  form  to  yield  a finite  result 
in  the  limit  It  -►  ®.  The  final  predictive  equation  is: 


a = 39.76  - 10.5  M + 91.6^4. 
m t 


standard  error  of  = 6.95 

m 


(15) 


Only  am  corresponding  to  Cym  > *75  are  represented  in  this  equation 
for  the  previously  mentioned  reasons,  and  for  data  having  two  Cym 
peaks,  both  are  used. 


The  variable  a©  is  related  to  It  by  both  the  data  base  and  the 
crossflow  model.  Accordingly,  the  following  predictive  equation  is 
pr-jposed: 


a 


o 


which  yields: 


(16) 


a = 18.22  + 

O -T  ^ 


standard  error  of  a =4.75 

o 


COMPARISON  OF  REGRESSION  AND  CROSSFLOW  MODEL 

The  crossflow  model  for  Cym  differs  significantly  from  the  final 
regression  equations . Although  both  appx'oaches 

indicate  that  M is  an  important  variable,  the  critical  model  dimension 
is  In,  not  It  as  suggested  by  the  crossflow  approach.  Also,  on  sharp 
tangent  ogives  peak  Cy  values  are  three  to  four  times  larger  than  the 
crossflow  model  predictions  and  only  a weak  Reynolds  number  effect 
is  noted.  These  differences  suggest  that  Cym  values  are  a product  of 
the  flow  about  the  model  nose.  The  expanding  body  radius  in  the  nose 
region  produces  a favorable  pressure  gradient  retarding  the  formation 
of  a turbulent  boundary  layer  and  reducing  the  Re  effect.  To  correctly 
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model  Cym/  Cl  values  should  be  used  which  correspond  to  the  initial 
ones  following  the  start  of  impulsive  flow  about  an  expanding  cylinder. 
Such  experimental  information  is  currently  not  available. 

The  prediction  of  the  crossflow  model  for  Cq  closely  corre- 
sponds to  the  results  of  the  linear  regression  model  as_is  shown  in 
Figure  14.  Both  approaches  define  o.^  as  a function  of  1^  which 
decreases  with  increasing  value  of  this  variable.  This  good  agree- 
ment indicates  that  the  onset  of  side  forces  occurs  in  a manner  which 
is  consistent  with  the  crossflow  model. 

The  crossflow  model  and  the  regression  equations  both  indicate 
that  ttju  is  a function  of  It  and  M/  decreasing  as  either  of  these 
variables  increases.  The  importance  of  Tt  in  the  regression  model 
is  unexpected  since  Ij^  is  the  critical  model  dimension  in  determin- 
ing Cyjn.  This  can  be  understood  by  noting  that  otuj  > oq  and 
Oq  varies  extensively  with  It-  The  physical  significance  of  the 
Mach  number  dependency  for  is  clearly  explained  by  the  crossflow 
model  which  indicates  that  Cl  and  hence  Cy  decreases  as  Me 
increases  towards  the  transonic  regime.  Hence,  for  M near  or  greater 
than  unity,  Cym  must  occur  at  low  a to  insure  that  M^  is  small. 

The  solution  for  ap  as  defined  by  Equation  (8)  is  plotted  in 
Figure  12.  This  quantity  represents  the  minimum  angle  of  attack  at 
which  a complete  cell  is  present  on  body.  As  can  be  seen  from  this 
graph,  it  seems  to  define  a lower  bound  for  am* 

CONCLUSIONS 

A detailed  examination  of  the  data  base  leads  to  the  following 
conclusions; 

(a)  Cones,  sharp  tangent  ogives,  and  blunt  tangent  ogives  all 
show  different  side  force  characteristics. 

(b)  In  all  cases  Cym  decreases  with  increasing  M and  disappears 
for  M > 1.4  on  tangent  ogives  and  M > 2.0  on  cones. 

(c)  On  tangent  ogives  Cym  decreases  with  increasing  In  partic- 
ularly in  the  case  of  sharp  bodies. 

(d)  Blunt  tangent  ogives  show  a lower  level  of  Cym  than  sharp 
ones,  but  the  degree  of  bluntness  does  not  appear  important. 

(e)  The  onset  angle,  oq,  decreases  with  increasing  It. 

(f)  The  parameter  Om  decreases  with  increasing  M and  it. 

The  final  regression  equations  incorporate  the  above  observations 
and  in  addition  suggest  the  following  which  is  motivated  by  the  cross- 
flow  model  predictions. 
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(g)  On  tangent  ogives,  Cym  is  a weak  function  of  Re,  decreasing 
under  the  conditions  of  turbulent  separation  at  subcritical  Mach 
numbers. 

This  last  conclusion  should  be  considered. as  tentative.  Most  of 
the  high  Re  points  are  from  a single  facility-’-and  the  supposed  Re 
influence  may  possibly  be  a facility  effect. 

Comparison  of  the  linear  regression  equations  with  the  crossflow 
model  predictions  suggests  that  the  peak  loads  on  sharp  bodies  are 
associated  with  tlie  flow  about  the  body  nose.  This  conclusion  is 
based  on  the  observation  that  is  the  critical  dimension  in  the 
regression  model  and  plays  a role  similar  to  1^  in  the  cross" low 
model.  The  good  agreement  between  regression  and  crossClw«  models 
for  olq  indicate  that  side  forces  arise  in  a manner  which  is  consistent 
with  the  crossflow  theory. 

The  linear  regression  equations  for  Cy^i,  Uq,  and  developed 
in  this  report  provide  for  the  first  time  a method  of  estimating 
values  for  these  parameters.  These  specific  equations,  which  should 
be  updated  when  a more  comprehensive  data  base  becomes  available,  are 
a product  and  not  a part  of  the  papers  theme.  The  central  idea 
is  that  realistic  prediction  methods  for  bodies  at  high  angles  of 
attack  in  the  subsonic -transonic  range  must  take  into  account  the 
nondeterministic  nature  of  the  flowfield.  The  load  on  and  hence,  the 
flowfield  about  bodies  under  these  conditions  will  fluctuate  for 
reasons  which  are  difficult  to  identify  and  often  not  practical  to 
control.  Conclusions  concerning  the  characteristics  of  side  force 
cannot  be  drawn  from  a limited  number  of  tests,  but  can  only  be 
determined  in  a statistical  sense  by  considering  a large  number  of 
measurements.  However,  a physical  understanding  of  the  phenomena  is 
essential  for  the  development  of  an  optimum  predictive  model. 
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